Pion-Photon Transition Form Factor and New Physics in the r Sector 



David McKeen/'B Maxim Pospelov, 1 - 2 'Q and J. Michael Roney 1 '!! 



(N 

o 

(N 



Victoria, BC V8P 5C2, Canada 
ON N2J 2W9, Canada 



1 Department of Physics and Astronomy, University of Victoria, 
2 Perimeter Institute for Theoretical Physics, Waterloo, 

(Dated: March 15, 2012) 

Recent measurement of the 77* form factor of the neutral pion in the high Q 2 region disagrees 
with a priori predictions of QCD-based calculations. We comment on existing explanations, and 
analyze a possibility that this discrepancy is not due to poorly understood QCD effects, but is a 
result of some new physics beyond the standard model (SM). We show that such physics would 
necessarily involve a new neutral light state with mass close to the mass of it , and with stronger 
than 7T° couplings to heavier SM flavors such as c, r, and b. It is found that only the coupling 
to the t lepton can survive the existing constraints and lead to the observed rise of the pion form 
factor relative to Q~ 2 at high Q 2 . We perform numerical fits to data and determine the allowed 
range of masses and couplings for such new particles. This range of masses and couplings could 
also reduce or eliminate the tension between the e + e~ and r decay determinations of the hadronic 
vacuum polarization. Dedicated experimental analysis of r pair production in association with such 
new states should provide a conclusive test of the new physics hypothesis as an explanation of the 
pion form factor rise. We also comment on the calculations of the pion form factor in the chiral 
quark model, and point out a possible dynamical origin of the quark mass scale inferred from the 
form factor measurement. 



PACS numbers: 13.40.Gp, 14.40.Be, 14.60.Fg, 14.80.Ec 



J2' 

> 

O 

(N 
(N 

(N 



X 



I. INTRODUCTION 

The studies of properties of light pseudoscalar mesons 
via the "photon fusion" mechanism has been a classic 
subject in particle physics that predates quantum chro- 
modynamics (QCD). With the advent of QCD came the 
impetus to study the off-shell form factors of 7r°, 77, and 77' 
as a way of testing QCD in exclusive processes. The most 
recent experimental advance in this direction has been 
achieved by the BABAR collaboration 0, Q that measured 
the form factors of pseudoscalar mesons in the region 
of momentum transfer extending to \Q 2 \ < 100 GeV 2 , 
significantly increasing the experimental Q 2 range over 
previous measurements by CLEO Q and CELLO [|J. 
The pion transition form factor measurement is based on 
the process e + e~ — > e + e~77* — > e + e~ir°, with one lep- 
ton scattered at a small angle and emitting a quasireal 
photon, and the other lepton scattering through a large 
angle, thus emitting a deeply virtual 7*. No sign of 
convergence towards perturbative QCD asymptotics [f| 
is seen in the BABAR data for the it while the r\ 
and r( form factors Q appear to be consistent with 
prior theory expectations. The probed momentum trans- 
fers, Q 2 ~ 40 GeV 2 , extend far beyond any nonper- 
turbative energy scale relevant for low-energy QCD - 
Q 2 > to 2 , Agj CD , m* W) , (47r/ w ) 2 , which is why the 
disagreement of experimental results for ir° with QCD- 
based theory expectations came as a total surprise. 

An intense theoretical debate of the significance and 
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interpretation of BABAR result ensued, see, e.g., Refs. 6 
Hoj , and on the basis of this discussion the following log- 
ical possibilities can be identified: 

1. The measurement of the tt° form factor by BABAR 
has some unidentified error. 

2. Pre-BABAR QCD predictions were somehow 
"naive," and need to be augmented by assuming, 
e.g., considerably different patterns of the quark 



momentum distribution within the 



n 



4. 



Alternative ways of describing the tt° transition 
form factor, such as, e.g., a constituent quark 
model, should be preferred over the QCD descrip- 
tion. 

The discrepancy of the BABAR measurement with 
QCD predictions is in fact a signature of new 
physics beyond the standard model (SM). 
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While first three options on this list have been widely 
discussed in the literature, the last possibility, new 
physics, looks the least likely and thus far has been ig- 
nored. Indeed, a new physics option would have to be 
rather contrived: it must involve new particles/states be- 
low the weak scale, have interactions with photons and 
7T° far stronger than the weak interactions, and evade all 
other constraints and tests. 

In this paper, we take this last option as an assump- 
tion and investigate broad classes of new physics capa- 
ble of influencing the BABAR. form factor measurement. 
We identify two classes of models that lead to signifi- 
cant modifications of the form factor. In the first model, 
a pion mixing with a new light state results in the en- 
hanced coupling of 7T° to intermediate mass quarks, c, b 
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and to the r lepton. We call this model a "hardcore" 
pion. In the second class of models a new particle, scalar 
or pseudoscalar, is introduced with its mass very close to 
the 7T° mass. Such new particles can be called "pion im- 
postors." We show that the loop contributions of heavy 
quarks and r leptons are capable of changing the form 
factor behavior, but only the coupling to r has a reason- 
able chance to remain undetected, while couplings to b 
and c quarks directly to tt° or to a pion-like impostor par- 
ticle are impossible to reconcile with bounds from bb and 
cc physics. We then identify important search channels in 
r pair production and decay which can be experimentally 
distinguished from SM backgrounds. Therefore, these 
models can be decisively tested with the future analyses 
of the r lepton dataset. 

This paper is organized as follows. Section |TT] intro- 
duces the basic formalism of form factor physics, dis- 
cusses the benchmark QCD predictions, as well as alter- 
native approaches within the chiral quark models. We 
choose the pre-BABAR QCD prediction [ll| that we use 
together with new physics contributions in our numeri- 
cal fits. Section ITO1 introduces models of new physics and 
includes calculations of corrections to the form factor. In 
Sec. IIV1 we determine the optimal ranges of parameters 
of the model, that is the strength of couplings to r lep- 
tons and the mass of the impostor particle. Section [V] 
describes further constraints on the models and possible 
strategies in detecting a new (pseudo) scalar particle pro- 
duced in association with r pairs and in r decays. Our 
conclusions are reached in Sec. [VI] 

II. PION FORM FACTOR IN QCD AND IN 
QUARK MODELS 

We begin by introducing the matrix element for the 
relevant reaction, 7*7* — > P, where P stands for a pseu- 
doscalar particle. It can be written as 

M = ie 2 T (q 2 , Q 2 ) e^e^k^ (1) 

with ei and £2 the polarizations of the photons, k\ and 
ki their momenta, k\ = —q 2 , k\ = —Q 2 , and J- (g 2 , Q 2 ) 
the transition form factor. Calculations and measure- 
ments of T have a long history, and the form factor for 
the 7r°-photon transition Q at large photon virtualities 
is among the earliest exclusive calculations in perturba- 
tive QCD. At large momentum transfers, the asymptotic 
freedom of QCD allows the matrix element to be factor- 
ized into the product of a hard scattering amplitude, T, 
calculated perturbatively, and a universal pion distribu- 
tion amplitude (DA), which encodes nonperturbative 
effects: 

F{q\Q 2 )=l dx T (q 2 ,Q 2 , fi 2 F , x) 4>n (x, fi%) (2) 
Jo 

where fip is the factorization scale and x is the momen- 
tum fraction carried by the quark in the pion. The DA is 



conveniently expanded in terms of Gegenbauer polynomi- 
als and obeys the Efremov-Radyushkin-Brodsky-Lepage 
(ERBL) evolution equation, with the asymptotic form 
4>^ sy (x,n 2 F ) — (v^/tt/S) 6x (1 — x) where is the pion 
decay constant. (We use the notation f n = V^F^ = 
130.4 MeV.) 

Experimentally, it is easier to access the region where 
one of the photons in the reaction is nearly on-shell, 
q 2 — > 0, rather than both photons having large virtuali- 
ties. Consequently, the large distance scale introduced by 
the on-shell photon necessitates use of nonperturbative 
techniques such as light-cone sum rules (LCSR) where 
F{q 2 iQ 2 } is related to a dispersion relation in q 2 , re- 
sulting in an expression for the form factor T (Q 2 ) 
.7 ' [cj 1 = 0, Q 2 ) in terms of an integral over a spectral den- 
sity. Despite this complication, the pion form factor at 
sufficiently large Q 2 should approach the Brodsky-Lepage 
limit, F BL = V2fn/Q 2 ^ 185 MeV/Q 2 . However, how 
exactly T reaches this limit as a function of Q 2 is open 
for debate, and many alternative ways of parameterizing 
the pion DA that could delay the onset of the asymptotic 
regime have been proposed over the years [Til fl2j . 

A reasonable agreement of CLEO [|| and CELLO [3] 
results with the QCD sum rule-derived pion form factor 
gave support to the overall consistency of QCD-based 
theory and experiment. A recent result from BABAR on 
the 77 and 77' transition form factors is also in agree- 
ment with QCD 0. In light of this, the BABAR result 
on the pion form factor increasing well above Brodsky- 
Lepage asymptotics and showing no sign of Q 2 F(Q 2 ) 
"flattening," definitely looks puzzling. This increase in 
the form factor has spurred a large amount of theoreti- 
cal work revisiting the issue. Notable post-BABAR QCD- 
based explanations one way or another employ a "flat" 
DA 0]. Also, serious doubts were voiced as to whether 
the BABAR-observed growth of Q 2 F(Q 2 ) is actually possi- 
ble within consistent theories of strong interactions [l(| ■ 
In our paper, we will not side with any of these claims, 
simply noting that none of QCD-based papers could pre- 
dict the pion form factor behavior, as observed by BABAR 
above Q 2 = 10 GeV 2 . In light of this, in our subsequent 
new physics speculations, we are going to use the QCD 
form factor as predicted by the most elaborate pre- BABAR 
QCD-based study in Ref. [HI- 

The stark disagreement of the measured pion form fac- 
tor with QCD predictions has motivated some authors 
to evaluate the same quantity within various quark mod- 
els. These models operate with "elementary" pions and 
quarks, coexisting as dynamical degrees of freedom. Con- 
sider a quark model where the derivative coupling of tt° 
to light quarks is governed by g a — 1 coupling, 

1 <9 t 7r° 
£qm = -d fl TT a d ll 'K a + qiiD^f, - m q )q + -^-^qinlh^q, 

which is a subset of the full quark- meson Lagrangian [L3| , 
and has explicit chiral symmetry in the limit m q — > 0. 
Typically, these models assume that m q is given by the 



3 



constituent quark mass scale, and therefore is not small. 
A simple calculation of a quark triangle diagram in the 
regime Q 2 3> rn 2 with two quark mass insertions gives 
well-known double-log asymptotics for the pion form fac- 
tor, 
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(4) 



The current quark mass values of O(10 MeV) provide 
a negligible contribution to the form factor at large Q 2 , 
but constituent quark masses on the order 0(100 MeV) 
can provide a reasonable fit to the data, as was pointed 
out in a number of publications [&, 8] . The naive asymp- 
totic formula of Eq. ((U) can be modified to account for 
the fact that beyond a certain nonperturbative scale A, 
the constituent mass m q is expected to ge t an additional 



power-like suppression, m q ~ A /Q 



and that the 



pion-quark vertex may also be cutoff at certain quark 
virtualities. All this will lead to a stabilization of the 
logarithmic growth, and in light of this the following phe- 
nomenological formula is preferred: 



T^a (Q 2 



1 



m 2 / log 2 (0 2 /m 2 ) for Q 2 < A 2 
Q 2 \ log 2 (A 2 /™ 2 ,) forO 2 > A 2 



( 5 ) 

The functional form of the transition between these two 
regimes is of course impossible to track without extra 
assumptions about the m q (Q 2 ) behavior. 

A fit to the experimental data can provide information 
on the preferred values for the parameters m q and A. 
To interpolate between the two regions of Eq. [5J we take 
Q 2 A 2 / m 2 (Q 2 + A 2 ) as the argument of the logarithm. We 
plot the results of such a fit in Fig.Q]where 3er contours of 
preferred regions are shown when all experimental pion 
form factor data are used (black, solid), and when only 
the BABAR results are employed (red, dashed). This fig- 
ure shows that when all data are included in the fit, the 
preference is for m q in the range from 130 to 200 MeV, 
and A between 2 and 5 GeV. When only BABAR data are 
included, there is a clear preference for larger A, includ- 
ing values approaching infinity, in which case the quark 
masses should be in a relatively narrow range around 
m q ~ 130 MeV. This is consistent with results found by 
other authors @. 

How legitimate are such explanations? Some authors 
Q remark that a mass scale of 130 MeV looks very odd: 
it is far above the current quark mass values, and yet 
significantly below a "typical" constituent mass value of 
^300 MeV. Moreover, if both light quarks and pions are 
included within the same dynamical scheme, one should 
be worried about the self-consistency of the theory at the 
loop level, renormalization of F„, etc. 

Here we would like to point out how a mass scale 
m q ~ 130 MeV could emerge in the theory described 
by ([3]). We consider the leading order power correction 
to the propagator of a light quark in QCD and in the 
chiral quark model. For simplicity, we take the value of 
the current quark mass m = (m u + m ( i)/2 close to zero 
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> 




FIG. 1. 3<r contours resulting from the fit of Eq. ((5| to the 
experimental extraction of pion form factor. The solid (black) 
contour corresponds to using all experimental data, and the 
dashed (red) to using only BABAR data. The best-fit point 



marked inside the black contour is at m q 
A ~ 2.6 GeV. 



160 MeV and 



(m w — » limit), while maintaining a finite value of the 
quark vacuum condensate (qq) ~ —(270 MeV) 3 . For the 
case of two light flavors with SU(2) symmetry, direct cal- 
culation show that the quark propagators can be written 



<Sqcd(p) 



i i4na s (qq) 



3p 4 



Scqm{p) = ]>- v&fIp 2 



(6) 
(7) 



where ellipses stand for higher-order power corrections in 
terms of the ratio of nonperturbative scale of strong in- 
teraction to p 2 . The power corrections can be matched to 
the S(p) — i/p' + im q s /p 2 , and the second terms in these 
formulae can be identified with the quark mass contribu- 
tion. The QCD result © is of course a well-known for- 
mula due to Politzer in the Landau gauge [13, HH (with 
the numerical error of [Tij corrected), resulting from the 
quark-gluon propagator mixing, Fig. [5] The second line, 
Eq. ([7]), is the result of the quark- meson propagator mix- 
ing in the chiral quark model. The scale of the effective 
quark mass suggested by (J7J is 



,eff 



(qq) = n 

16F 2 ~ 16(m„ 



135 MeV x 



4.3 MeV 



in 



where the current quark mass value is within the range 
of 3.0 — 4.8 MeV at a normalization poin t of 2 GeV, as 
quoted by the Particle Data Group [16j. Interestingly 



enough, the mass scale emerging from the quark conden- 
sate correction to the quark propagator (J8J) is exactly in 
the range suggested by our fit to the pion form factor. 
If instead one considers three flavors of light quarks in 
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the chiral limit with flavor SU(3) symmetry, the result 
changes to (qq)/(9F%) ~ 240 MeV. We stress that © is 
an effective quark mass that neither breaks chiral sym- 
metry explicitly nor makes the pion massive. 

The QCD result is clearly inapplicable below 1 GeV, 
and moreover, it has an explicit gauge dependence, that 
will cancel in the full one-loop evaluation of the gauge- 
independent correlators of quark currents [15I | . The result 
for S'cqm is °f course gauge independent. Because of 
this distinction, the matching of S'cqm and Sqcd cannot 
be done rigorously. Nevertheless, the power corrections 
in © and ([7]) become comparable around p 2 ~ 1 GeV 2 , 
which is precisely the scale where one would expect the 
quark model description to become invalid and be re- 
placed by perturbative QCD. One can see that this scale 
is a lot lower than scale A suggested by the numerical fits 
to the form factor. 

The closeness of in Eq. ([8]) to the value suggested 
by the chiral quark model interpretations of the pion 
form factor is intriguing. However, at this point it cannot 
be considered as a legitimate explanation of the EABAR- 
observed rise of the pion form factor. One has to use the 
chiral quark model well beyond \p 2 \ ~ 1 GeV 2 and indeed 
outside of its usual range of validity. Moreover, even as- 
suming that quark-pion Lagrangian can be extended that 
far in energy, one would still need to calculate the chiral 
loop corrections to the triangle diagram to see at what 
energy scale the whole framework becomes inconsistent. 
Within SU(3) flavor symmetry, it is also difficult, if not 
impossible, to explain a markedly different behavior of 
the form factors for the pion and the SU(3) octet com- 
bination of r\ and rj [Hi . Nevertheless, the emergence of 
the effective quark mass due to the quark condensate, 
Eq. ((5J, within massless chiral quark models calls for 
the evaluation of various quark current correlators in the 
same framework. A procedure very similar to the QCD 
operator product expansion (OPE) can be developed this 
way, with the hope that ([3]) could provide a reasonable 
description of hadronic correlators in the intermediate 
range of energies, (m° ff ) 2 <C |j> 2 | <C GeV 2 . We will leave 
the investigation of the hadronic models with pions on 
the "wrong side" of the OPE to future work. 

To summarize this section, QCD calculations of the 
pion form factor were remarkably unsuccessful in predict- 
ing the rise of the pion form factor observed by BABAR. 
Chiral quark models may provide a reasonable explana- 
tion, but have to be extended far beyond the usual limits 
of their validity. 



III. NEW PHYSICS 

In this section, we describe an explanation of the rise 
in the 7r°-7 transition form factor using new physics. We 
will consider r leptons coupled to either new pseudoscalar 
or scalar degrees of freedom with masses close to that of 
the neutral pion. In the case of a new pseudoscalar, there 
is also the possibility of mixing with the ir° so that an 




FIG. 2. Nonperturbative corrections to the quark propagator 
in QCD and in the chiral quark models. 

effective coupling of tt° to r is induced. As will become 
clear later, new neutral particles coupled to r-lcptons is 
perhaps the only model of new physics that can provide 
an appreciable contribution to the pion transition form 
factor and evade other constraints. 

The masses of the new light states we consider are es- 
sentially free parameters in this phenomenological study 
and only fixed according to experimental constraints. 
Their small masses do not necessarily result from their 
being psuedo-Nambu-Goldstone bosons of a sponta- 
neously broken global symmetry as in the case of the 
pions. A precise embedding of the models we consider in 
a more complete theory could, however, offer an expla- 
nation. 



A. A New Pseudoscalar 



We first consider a SM singlet pseudoscalar, P, coupled 
to SM fcrmions through dimension-5 operators, 




(9) 

where H° is the neutral component of the SM Higgs dou- 
blet and A is the scale at which this effective interaction is 
generated. The ellipsis refers to other couplings, e.g., to 
other SM fermions, which we will ignore for the time be- 
ing. Equation © is the unitary gauge expression for the 
SUl(2) x UY(l)-invariant effective interaction. When the 
Higgs field is replaced by its vacuum expectation value 
v/y/2, effective Yukawa interactions result, 

Ant = -iP (<?tT7 5 t + g u u-/ 5 u + g d dj 5 d) , (10) 

with g u = \ u v/\/2A, etc. Since we expect A to be near 
the scale of electroweak symmetry breaking, the effec- 
tive Yukawa couplings g u ,d,T might not be small unless 
required by experiment. 

An axial-vector coupling and a photon coupling 
are available for the pseudoscalar particles as well: 
/^ 1 <9 A1 Pf7 M 7 5 T and f^PF^F^. For each flavor of 
fermions these three couplings are related by one anomaly 
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relation. We choose to eliminate the axial-vector cou- 
plings and keep the pseudoscalar and photon couplings. 
Moreover, the "contact" photon coupling is not suit- 
able for fitting the pion form factor, as it gives a Q 2 - 
independent contribution to F{Q 2 ) to be compared with 
the fiducial 1/Q 2 fall off. 

We can now imagine two scenarios distinguished by 
whether there is mixing between P and tt°. 



We will refer to P in this scenario as a "pseudoscalar im- 
postor." In this formula, we also assume that the contact 
photon coupling of P to F^F^ is small. If it is not, then 
the Fp(Q 2 ) contribution in (fTBj) will receive a constant 
Q 2 -independent addition. 

2. A Hardcore Pion 



1. Pseudoscalar Impostor 

First, if the Yukawa couplings of P to u and d are 
vanishing, mixing between P and 7r° is negligible, and 
the only interaction involving P and SM fields that is 
important is 



Ant = -igrPr^r. 

This leads to the following amplitude for 77* 
through a r loop, 



M. p 77 * —ie 2 Fp 



fc l t 2"'l ™2 > 



(11) 
-> P 

(12) 



where k\ and ki are the photons' momenta with k\ — 
—q 2 = and k\ = —Q 2 . The form factor here is given 
by a well-known result for the fermion triangle diagram 
181 



Fp (Q 2 ) = 



g T m T 



7T 2 (Q 2 + m 2 p ) 



-1 rn P 



2to 7 



with 



L = — log 

2 B 



R02 - 1 



R, 



Q 2 



1 



4m 2 . 

~Q 2 ~ 



(13) 



(14) 



If kinematically allowed, P would be produced in 



77 



z + e~P and primarily decay to 
77. Further, if the mass of the pseudoscalar is tuned 
so that mp ~ m^o (a quantitative study of the extent to 
which this must hold will be given in Sec. IIVB[) then P 
would be misidentified as a ir° when reconstructed using 
a peak in the M 77 spectrum. The ir a -j form factor is 
extracted from the measurement of the cross section for 



e + e 7T° using 



da (e + t 



dQ 2 



\F(Q 2 )\ 
Q 2 



-A(Q 2 ) (15) 



where A (Q 2 ) is a function that depends on the particu- 
lars of the experiment and on radiative corrections. The 
matrix element in Eq. f| 1 2[) leads to a cross section for 



e+e' -> e+e-P given by 7ra 4 (|j> (Q 2 )| 2 /Q 2 ) A (Q 2 ) 

which, since mp ~ m^o, is interpreted as an additional 
contribution to the 7r° production cross section. This 
would lead to a modification of the form factor extracted, 



\F(Q 2 )\ 



F(Q 2 )\ 2 + \F P (Q 2 )f 



(16) 



The second scenario that could result from the inter- 
action in Eq. (|10j) is if the Yukawa couplings to u and 
d quarks are large enough to lead to appreciable mixing 
between P and 7r° if, again, mp~ m^o . 

In this case, an effective interaction of the 7r° itself with 
the r is generated, 



A 



-ig^f^r, 



(17) 



where <j£ = g T sin# ~ g T 9, and the effective mixing angle 
between P and 7r° is given by 



I (««) I (ffu - 9d) (400 MeV) 2 x (g u - g d ) 



i^(m 2 



i 2 P ) 



2m„.o AmpjO 



(18) 



We shall take the mass splitting to be small, AmpjO = 
Im^o — mp| <C mjr. 

As in the no-mixing case above, besides a correction to 
77* — > 7T°, there will also be an amplitude for 77* — > P, 
potentially complicating the analysis. To keep things 
simple, in this case with mixing between P and 7r°, we 
can make the production of real P's experimentally unim- 
portant by adding a long lived singlet state, <fi coupled to 
P, so that P decays invisibly, P —> tficf). To avoid making 
the invisible width of the ir° too large, we can forbid the 
decay 7r° — > <jxf> kinematically if, e.g., the mass of 4> is 
70 MeV and mp = 145 MeV. This way, the only change 
in the apparent form factor in the mixing case is given by 
Eq. (Ti"9|) and we do not deal with the added complication 
of a second state in the 77 spectrum just as we do in the 
no-mixing (pseudoscalar impostor) case. 

Given the interaction in Eq. (I17[) . there will be a change 
in the actual ir° form factor. This change is simply 



J" (Q 2 ) -> T (Q 2 ) + F P (Q 2 



(19) 



where Fp (Q 2 ) is given by Eq. (|13l) with the replacements 
g T — > mp — > m^o . We will call this the "hardcore 
pion" scenario. 

Because the hardcore pion model uses couplings of new 
states to light quarks it should be considered "more ex- 
treme" than the pseudoscalar impostor model, and is 
subject to additional constraints coming from the quark 
sector. Therefore, it deserves special discussion about its 
viability. The angle 9 characterizes the admixture of the 
new P state in the physical ir°. As such, the well-tested 
rate for the decay 7r ± — > 7r°e ± z/ will be corrected, 



r(7r 



) = r SM x cos 2 9 ~ r SM x (l-9 2 /2). 

(20) 
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Prom the 1% agreement between theory and experiment 
for this rate one should conclude that 9 2 < 0.02. To sat- 
isfy this limit with mp = 145 MeV, the Yukawa couplings 
should be at the level of \g u — gd\ < 2 x 10~ 3 , which is not 
a very tight bound. The mass shift of the physical pion 
is ~ Amp^oO 2 , and is less than 0.2 MeV if Amp^a < 10 
MeV, within the uncertainty of the theoretical calcula- 
tion of the neutral-charged pion mass difference (see, e.g., 



Since we do not impose any symmetry requirements 
on g u .d, then one should also expect additional viola- 
tions of isospin invariance due interactions mediated by 
P. The measure of such new effects is on the order of 
9u d/( m ?/^) 2 ' which again should not exceed ~ 1%. 
One can see that with g u ^ ~ O(10 -3 ) such constraints 
will be easily satisfied. The effects of a possible mixing 
with r\ are not important, as we assume that P is much 
closer in mass to 7r°. 

Another important constraint could come from missing 
energy decays of K mesons, K + — > ir + P —> n + f!, since 
we assume that P decays invisibly. The branching ratio 
for this decay is tested to O(10~ 10 ) [2(J, and indeed could 
provide a powerful constraint, except in the region of 
parameter space with \m^o — mp\ < 10 MeV, where the 
measurement of the missing energy rate is not possible 
due to 7T + 7r° background. 

We conclude that based on the arguments of isospin 
invariance, precision tests of charged pion beta decay, 
and the missing energy decay of K + , the mixing angle 9 
is limited to ~ 0.15, while the mass of a new pseudoscalar 
state should be close to the pion mass within ~ 10 MeV. 
Finally, the effect of this mixing on the rate for ir° — > 77 
is discussed in Sec. IV Bl 



B. A New Scalar 



If we consider a new scalar instead of a pseudoscalar, 
the analysis is similar but simplified. Since the scalar 
mixing with tt° breaks CP symmetry, we ignore couplings 
to light quarks and only consider couplings of the scalar 
to r through 



Ant = — & ( 



(A' r fr) -> -hrSfT. 



(21) 



As in the pseudoscalar case, this generates an amplitude 
for two photons to produce a scalar, 77* — > S, 



M Sll =e 2 T s (Q 2 ) {k»k% - h ■ k 2 g fiu ) e lM e 2 „ (22) 



and we have again taken k\ = —q 2 = 0. The form factor 
in this case is 



h T m T 



1 - 



tt 2 (Q 2 + m 2 ) 
Am 2 



Q 2 



L 



-1 m s 



2Q 2 



Rq2 L + R m i ( sin 



with L and Rq2 as in Eq. ([T4|) and 



2m T 

-1 ms 
2m T 



(23) 



/ Am 2 



(24) 



Defined in this way, the e + e — > e + e S cross section is 



da (e + e —> e + e S) 
dQ 2 



(Q 2 
Q 2 



-A (Q 2 ) (25) 



with A (Q 2 ) the same as in the pseudoscalar case. There- 
fore, if the scalar mass is close to that of the neutral pion, 
it could mimic a 7r° giving an apparent contribution to 
the 7r°-7 form factor, 



\T(Q 2 )\^^/\T(Q 2 )\ 2 + \T S (Q 2 )\ 2 . (26) 
This particle will be referred to as a "scalar impostor." 

C. Why r? 

We now turn our attention to the motivation for choos- 
ing r in the interactions in Eqs. (jTTj) . (TT7|) . and (|21l) . 
First, if we wish to have a change in the (apparent) 
7T°-7 form factor due to a loop-induced coupling of an 
impostor or the pion to 77, the new states running in 
the loop must obviously be electrically charged. If it is 
a state unaccounted for in the SM, then its mass must 
be at least on the order of 100 GeV to escape collider 
bounds on pair production of new charged particles. Par- 
ticles with masses this large will require extremely large 
couplings to P or S to have any effect on F (Q 2 ) for 
Q 2 ~ 10 — 40 GeV 2 . Thus, we are led to only consider 
charged SM fermions (excluding the top quark) as can- 
didates for coupling to impostors or to the nr. 

In order to change the form factor behavior, the cou- 
plings of new states to light quarks would have to be very 
large, stronger than the pion coupling to quarks, which 
is clearly unrealistic. Therefore, we ignore couplings to 
u,d,s, only considering small values that potentially me- 
diate mixing with the pion in the hardcore pion case but 
are otherwise negligible. 

Couplings to charm and bottom quarks could explain 
the rise in the pion form factor but run into limits from 
the decays of cc and bb states. 
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In the charm quark case, the couplings can be con- 
strained by data on J/tp — > 771"° and ip' — > jir° since 
pion impostors that decay to 77 will be reconstructed as 
neutral pions. Requiring that B (?/>' — > 77r°)not exceed 
its experimental value of (1.58 ± 0.42) x 10~ 6 [Il| sets 
limits on the coupling of impostors or a hardcore pion to 
c quarks of 

\g c \,\9c\ £ 0.018, \hc\< 0.023 (27) 

which are too small to have an appreciable effect on the 
7T° form factor. 

To constrain the coupling to bottom quarks, we note 
that b quarks interacting directly with the ir° or a scalar 
or pseudoscalar impostor will lead to (apparent) enhance- 
ments in the transitions between bb states with the emis- 
sion of a 7T°. Approximating the interquark potential 
as being Coulombic and using the experimental limit, 
B(T(2S) T(1S>°) < 1.8 x 10~ 4 at 90% C. L. [Hj], 
results in the limits on the couplings to b quarks of 

\g b \, \g%\ < 0.055, \h b \ < 0.0055. (28) 

As in the case of the charm quark, these are too small to 
give a non-negligible contribution to the form factor. 

Moving on to leptons, we can immediately rule out 
the electron as a candidate since this will lead to 77 no 
longer being the dominant decay mode for P or S in the 
impostor cases, and in the hardcore pion case leading to 
a branching ratio for ir° — > e + e~ that is many orders of 
magnitude larger than experimentally allowed. 

The strictest limits on couplings to muons come from 
(g — 2) M . These couplings must be suppressed to the 
level oflO -3 to 10 ~ 4 , and therefore coupling to the muon 
cannot offer an explanation of the rise in the form factor. 

For these reasons, we concentrate on the scenario 
where the ir° or an impostor is coupled to r. 



D. Limits from (g — 2) T 

The most immediate constraints on the interactions 
of the t with light spin-0 particles, as in the case of 
the muon, are due to limits on the its anomalous mag- 
netic moment but, since it is much less precisely known 
than that of the muon, the constraints are weaker. Pseu- 
doscalars or scalars whose masses are equal to m T o inter- 
acting with t as in Eq. (|11[) or (|21l) will induce contribu- 
tions to a T = (g- 2) T /2 of [H[ 

(Aa T ) p = -6.21 x 10- 3 gl, (29) 
(Aa T ) s = 1.66 x lO~ 2 h 2 T . (30) 

The hardcore pion scenario results in the same contri- 
bution as in (Aa T ) p with g T — > g^. DELPHI provides 
the most stringent experimental limits on a T : —0.052 < 
a T < 0.013 at 95% confidence level [H. Taking the SM 
expectation into account and requiring that the total a T 



not fall outside of this range gives the following limits on 
the couplings: 

\g T \,\g* T \ <2.9; \h T \ < 0.84. (31) 

These limits are mild enough that coupling scalars or 
pseudoscalars to r can offer an explanation for the rise 
in the 7r°-7 form factor above the asymptotic value in 
QCD. We will postpone discussion of other constraints 
on coupling spin-0 particles to r such as higher-order 
corrections to (g — 2) ei/i , limits from Z — > jir and Z — > 
7T + r~, and r decays until Sec. fVl 

IV. FITTING THE tt -7 FORM FACTOR DATA 

A. Couplings 

To explore whether the interactions in Eqs. (fTTj) . (|17l) . 
or (|2"Tj) can explain the rise of the 7r°-7 form factor in 
agreement with behavior seen by BABAR, we first have 
to choose a description of this form factor in QCD. For 
the purposes of our fits we take the form factor obtained 
by a calculation of the DA using QCD sum rules with 
nonlocal condensates by Bakulev, Mikhailov, and Stefa- 
nis (BMS) [ll[ . BMS obtained a range of allowable form 
factors by varying the cutoff of an integral of a pertur- 
bative spectral function by ±10% in their sum rule. We 
allow the BMS form factor that we use in our fits to vary 
in this range. 

In the case of the pseudoscalar impostor we fit the data 
on T (Q 2 ) from BABAR, CLEO, and CELLO using [see 
Eq. GSJ] 

\? (Q 2 ) I = (Q 2 )\ 2 + \Fp(Q 2 )\ 2 (32) 

with Jbms the BMS form factor described above, which 
accounts for the QCD contribution, and Tp from Eq. (|13[) 
with rap = m^o . J-p depends very weakly on mp chang- 
ing by less than 1% for m^o /2 < mp < Im^o in the rele- 
vant Q 2 range. For the scalar impostor we do the same, 
replacing Tp with Ts from Eq. (|23[) with ms — m^o . We 
fit [Eq. CED] 

\F (Q 2 ) I = l-^BMS (<3 2 ) + ^HC. (<3 2 ) I (33) 

for the hardcore pion, where J-hCtt is given by Eq. f| 13[) 
with g T — > gl and mp — m^o. 

The results of these fits, along with the range of BMS 
form factors used, are shown in Fig. El The range of BMS 
form factors corresponds to a deviation from the data of 
at least 3.5a. 

In the pseudoscalar impostor case we find a coupling 
g T = 0.63^Q 22- The best-fit coupling deviates from the 
data at 0.8cr level and the total form factor using this 
coupling is shown by the solid line in Fig. [3] 

The fit in the scalar impostor case gives h T — 0.97^o'34- 
Much of this region conflicts with the constraint from 
(g — 2) r ; taking this into account we arrive at a best-fit 
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FIG. 3. The best-fit total form factors in the pseudoscalar 
impostor (solid), scalar impostor (dashed), and hardcore pion 
(dotted) cases when fit to the data from BABAR [J], CLEO 
and CELLO @] using the form factors of Eqs. and (135)) . 
The resultant couplings are g T = 0.63, h T — 0.84, and = 
0.18, see Table [TJ The shaded region indicates the range of 
BMS form factors. The scalar impostor coupling is limited 
by the (g — 2) T constraint [24| . 



coupling h T = 0.84 which fits the data at the 0.9a level. 
The total form factor with this coupling is shown by the 
dashed line in Fig. [3J 

Fitting the data in the hardcore pion case results in 
g^ = O.lS^o!^- The best-fit coupling here, which corre- 
sponds to a 1.6cr deviation from the data, is smaller since 
there is constructive interference with the BMS form fac- 
tor which describes the data quite well at low Q 2 . In 
Fig. [3J we show the total form factor using the best-fit 
coupling in this case with the dotted line . 

The pseudoscalar impostor fits the data best of the 
three scenarios with the scalar impostor running into 
some tension with the limits from (g — 2) T and the hard- 
core pion being constrained by the good agreement of the 
BMS form factor with the low Q 2 data. 

The situation is not drastically different if we restrict 
our fit only to larger Q 2 data. For Q 2 > 8 GeV 2 , the 
range of BMS form factors deviates by at least 4.9<r from 



the data. Fits result in g T = 0.70 



+0.17 , 
-0.25' ""r 



1.08 



+0.26 
-0.38' 



and gl = 0.25_q;^4. Again, h T = 0.84 is chosen as 
the best-fit coupling in the scalar impostor case due to 
(g — 2) T constraint. The best-fit couplings deviate from 
the data at the 1.1a, 1.5c, and 1.2a levels for the pseu- 
doscalar impostor, scalar impostor, and hardcore pion 
cases, respectively. 

We summarize the best-fit couplings and the level of 
their agreement with the data in Table |TJ 



B. Impostor Mass 

It is clear that in order for pseudoscalar or scalar im- 
postors coupled to t to offer an explanation of the exper- 
imentally seen rise in the 7r° form factor as in Sec. IIV A[ 
they must promptly decay to 77. Also, m° — mp : s must 
be sufficiently small so that these new particles are re- 
constructed as neutral pions in the reactions e + e~ — > 
e + e~(P, S). This will create an apparent contribution 
to the cross section for e + e~ — > e + e~7r°, which is used 
to extract the form factor. We would like to understand 
just how small m\ — m P,S is required to be. 

Using the best-fit couplings in the impostor cases 
above, along with the expressions for the cross sections 
for e + e~ — > e + e~(P, S, 7r°) as functions of the form fac- 
tors for 77* — > (P,S,ir°), we can estimate the number 
of actual pions and impostors produced in a given Q 2 
range. For instance, using the BABAR experimental setup 
and the best-fit coupling for the pseudoscalar impostor, 
g T = 0.63, implies that the ratio of the numbers im- 
postors to neutral pions produced is Np/N^a = 5.4% 
for Q 2 — 4 — 4.5 GeV 2 . This ratio is quite insensitive 
to the mass of the pseudoscalar, varying by less than a 
factor of 10~ 2 as mp ranges from m^o/2 to 2m 7r o. We 
then fit the M 77 spectrum obtained by BABAR in this 
U range (Fig. 6(a) of Ref. Q), which is used to esti- 
mate the 7r° production cross section, with two Gaus- 
sians, one centered at m^o and the other at a value of 
mp that we choose, on top of a quadratic background. 
We fix the widths of both to 9 MeV, set by the experi- 
mental photon energy resolution. Their heights are then 
related by Np/N^o (5.4% here). We repeat these fits 
while scanning mp over a broad range around m T o and 
record the change in x 2 from the value obtained when 
mp = m„o. The allowed range of mp will be given by 
masses that do not lead to unacceptably large \ 2 when 
fitting the measured M 77 spectrum. We then perform 
the same fits of the BABAR data for Q 2 = 11 - 12 GeV 2 
and Q 2 = 20 - 25 GeV 2 where N P /N„o = 24%, 62% 
respectively using the spectra from Figs. 6(b) and (c) of 
Ref. [l|. The ratios N P /N n are similarly insensitive to 
the particular mp as for Q 2 = 4 — 4.5 GeV 2 . 

The results of these fits are shown in Fig. [U The 
Q 2 = 20 — 25 GeV 2 data is most constraining on the al- 
lowed range of mp although there are fewer events than 
in the lower Q 2 ranges, due to the increased ratio of pseu- 
doscalar impostors produced. There is a roughly ±5 MeV 
range near m^a where the x 2 is n °t increased by more 
than a few in each Q 2 range. We note that there are 
51 — 56 degrees of freedom in each fit. 

We repeat this procedure for the scalar with its best- 
fit coupling of h T = 0.84 and show the changes in x 2 
in Fig. [5J The ratios of impostors to pions we use is 
Ns/N^o = 4.2%, 18%, 46% for Q 2 = 4 - 4.5, 11 - 
12, 20 — 25 GeV 2 respectively. There are fewer impostors 
produced in this case compared to the pseudoscalar one, 
as can be seen in the smaller form factor in Fig. [31 and so 
the constraints on m$ are less tight with x n °t increas- 
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\cg U Vjrc V 


\u£ ^ O Vjrc V 


couplingila x^/d-o.f. agreement 


couplingila xVd-o.f. agreement 


p 

s 

HCP 


0.63^22 35.5/35 0.8a 
°- 97 io.34 °- 84 37.8/35 0.9a 
0.18^012 45.3/35 1.6a 


0.70io26 10.8/9 1.1a 
i-OSiass °- 84 13 - 4 /9 l-5a 
0.25±g;?| 11.9/9 1.2a 


BMS 


72.7/36 3.6a 


49.3/10 5.1a 



TABLE I. The ±la range of couplings (p T , /i T , from Eqs. |T7|), and (|21|l . in the pseudoscalar impostor (P), scalar 

impostor (S), and hardcore pion (HCP) cases when fit to the form factor data in the Q 2 ranges indicated. We also 

show the xVd.o.f. and the corresponding level of agreement with the data for the best-fit value in each case. Note that the 
best-fit value of coupling in the scalar impostor case is limited to the largest fit value (0.84) compatible with the (g — 2) T 
constraint (24l]-the x 2 reported here is with respect to this value. Also shown is the level of agreement between the data and 
the BMS form factors [111 ] for each Q 2 range. The total 7r° transition form factors using the best-fit coupling in each case are 
shown in Fig. [3] 




100 110 120 130 140 150 

m p (MeV) 

FIG. 4. Constraining the pseudoscalar impostor's mass: the 
change in \ 2 fitting the M 77 spectra of Fig. 6 of Ref. [l[ to 
Gaussians centered at m n o and nip with a quadratic back- 
ground, relative to the \ 2 f° r m P = m ^o. The solid curve is 
for Q 2 = 4 - 4.5 GeV 2 , the dashed for Q 2 = 11 - 12 GeV 2 , 
and the dotted for Q 2 = 20 - 25 GeV 2 . The heights of the 
two Gaussians used in the fits are related by the values of 
Np/N^o in each Q 2 bin (5.4%, 24%, and 62%). There are 
51 — 56 degrees of freedom in each fit. 



ing by more than a few over a range of about ±8 MeV 
around m^o . 



V. FURTHER CONSTRAINTS AND TESTS 

As seen in Sec. IIV1 explaining the deviation between 
the measured ir°-j form factor and most expectations in 
QCD with a coupling of the r to n° or to a new pseu- 
doscalar or scalar with a mass within about 5 MeV of m w o 
requires relatively large couplings, of the order 0.1—1. We 
would like to see if such interactions are constrained by 
different measurements and, if not, how to probe them. 




- 1 1 ., . . 

-20 r 1 

100 110 120 130 140 150 

m s (MeV) 

FIG. 5. Similar to Fig. [4] but in the scalar impostor case. 
N s /N w o = 4.2%, 18%, 46% for Q 2 = 4 - 4.5, 11 - 12, 20 - 
25 GeV 2 . 



A. Anomalous Magnetic Moments 

We have seen in Sees. IIIIDI and IIVI that the current 
experimental limit on (g — 2) T is not very constraining in 
the case of a pseudoscalar impostor or a coupling between 
r and 7T°. However, it does limit the scalar impostor 
coupling to r to h T = 0.84, causing it to not fit the rise 
in the ir°-j form factor as well as the pseudoscalar. 

In addition to the one-loop contribution to (g — 2) r , 
interactions like those in Eqs. (|11[) . (|17l) . and (|21[) would 
contribute to (g — 2) e . M at the three-loop level suppressed 
by at least a loop factor ~ 1/16-7T 2 times a coupling 
squared relative to the leading two-loop contribution in- 
volving a t loop. 

The leading r contribution to a M is roughly 4.2 x 10 -10 , 
comparable to the current experimental and theoretical 
uncertainties [25| . Therefore, the three- loop contribution 
is safely below the uncertainty and we conclude that 0(1) 
couplings of the r to a hardcore pion or impostors are not 
ruled out by corrections to (g — 2)^. 
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In the case of the a e , the leading r contribution is 
down by a factor ~ m^/m^ relative to the muon case. 
Thus, the three-loop contribution to a e from an impostor 
or hardcore pion coupling to r is well below the upper 
limit of 15 x 10 ~ 9 from other precision measurements of 
a remaining consistent with a e |26j |. 

In passing, we note that a ir° transition form factor 
that rises relative to Q~ 2 at the level implied by the 
BABAR data does not contribute to a meaningful en- 
hancement in the hadronic light-by-light contribution to 
(5-2) M [13. 



B. 7T° Width 

Currently, the most precise measurement of the 7r° — > 
77 rate, to 2.8%, has been made by the PrimEx exper- 
iment [28[ , using the measurement of the tt° production 
cross section through the Primakoff process 7 Nuc. — > 
7T°Nuc. which can be related to T (tt° —> 77). A hard- 
core pion coupling to r or impostors with masses close 
to that of the ir° could lead to an apparent deviation 
between this lifetime measurement from theoretical cal- 
culations. 

Using the limit for the QCD tt° form factor at Q 2 = 
from the chiral anomaly, and the expression for the 
hardcore pion form factor in Eq. [SI we find a fractional 
change in the rate for ir° — > 77 of 

V^*g^(l + £±-)*0.1g;. (34) 



C. Z Decays 

A hardcore pion or impostor coupling to r would also 
lead to (or mimic) Z — >• jir° decays through a r loop. 
However, because of the suppressed vector coupling of 
the Z to t, these processes would contribute to this decay 
mode with a branching less that 10~ 8 for couplings g T , 
h T , g™ equal to unity, well below the current limit of 

5.2 x 10~ 5 m. 

An associated emission of impostor particles or hard- 
core pions can be also constrained by the non-observation 
of Z — >• t + t~7 events by OPAL 31] with 7 well- 
separated from r leptons. Since at high energy the signal 
from pion-like particles will be indistinguishable from 7, 
we calculate the branching ratio of the Z to r pairs plus 
an extra particle, limit its energy to be Ep g ^o > 1 GeV, 
require that cos < 0.95 where is the smaller of the 
angles between (P, S, 7r°) and r + or t~ , the constraint 
B(Z — » t + t~(P, S,tt )) to be less than the experimental 
bound on the branching to t + t~j of 7.3 x 10~ 4 . The 
resulting sensitivity to coupling constants is 

\gr\M<1.0; |^| < 1-1- (38) 

In order to convert this to strict constraints we need to 
choose energy Ep g ^a, where ir° and 7 signals cannot 
be differentiated. Choosing this energy to be ~ 10 GeV 
further relaxes the limits in Eq. (f38|) by a factor of a 
few. We conclude that such limits are not quite sensitive 
enough to test the couplings required for a pion impostor 
or hardcore pion to explain the form factor measurement 
of BABAR. 



Using the best- fit value g™ = 0.18 gives a change of about 
1.8%. If we use the value of obtained from fitting only 
the Q 2 > 8 rmGeV 2 data it becomes 2.5%. These shifts 
are smaller than the experimental uncertainty as well as 
the uncertainty of theoretical calculations, and thus these 
values of g™ are not ruled out. 

In the case of a pseudoscalar impostor, the apparent 
shift in the lifetime is simply 



r 



P— >77 



1 



gW m l° 



*77 x 7r u — ^77 

~ 2.6 x 10" 



64:TT 3 m 2 



"gl 



(35) 
(36) 



In the case of a scalar impostor this is instead 



<$r„ 



+77 



7T U — V77 



r„ 



+77 



1447r 3 mJ 



1.2 x 10~ 3 ^. (37) 



Using the best-fit values for these couplings leads to shifts 
that are below the experimental uncertainty and we con- 
clude that they are also allowed. 

Lastly, following [29(, we do not expect a large contri- 
bution to the rate for ir° — >■ e + e~ in any of the scenarios 
we have considered. 



D. Searching for an Impostor 

1. t Decays 

One obvious place to look for a hardcore pion or pion 
impostors that interact with the r as in Eqs. (|lll) . (|17l) . 
and (|2ip is in r decays. It is well known that hadronic 
decays of r can be used a way of extracting the hadronic 
vacuum polarization, and that r-based extraction gives 
larger result than the direct e + e~ — > hadrons by about 
a few percent [32| . This mild tension provides additional 
motivation to search for models with additional "hadron- 
like" particles coupled to r-leptons. 

Our models will lead to unusual decay channels such 



as t- 



■K°e ± Dv or r ± -> P£ ± uu and t ± -> S^vv, the 



latter two mimicking the first. Here I denotes final state 
electrons and muons. The branching ratios are estimated 
using CalcHEP [33( and found to be 



) ~3.5x 10- 5 (^) 2 , 



B (t ± ->■ Pe ± vv) ~ 3.5 x 10~ 5 .g^ 



BlV* -> Se ± vv\ 



1.7 x IQ-^ht 



(39) 
(40) 
(41) 
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The best-fit values of the couplings lead to 



B (t ± Ti e ± Vv) ~ 1.1 x 10" 6 , 
B (r+ -> Pe ± vv) ~ 1.4 x 10~ 5 , 
B (r+ -> Se ± vv) ~ 1.2 x 1CT 3 . 



(42) 
(43) 
(44) 



The sum of all tau branching fractions, 0.9984 ± 
0.0010 [34], although consistent with unity does allow 
for B (t ± -> (P,5,7T°)^ ± w) < 3.6 x 10~ 3 at 95% confi- 
dence level. Such modes would not have been detected 
in conventional r branching fraction measurements and 
the related tau hadronic modes that include (P,S,ir°) 
would not have been detected as they would be included 
unknowingly in conventional hadronic branching fraction 
studies. Only a dedicated search for r+ — » (P, S, tt°)£ ± Du 
would reveal its presence in r lepton decay and to the 
best of our knowledge no such experiment has been per- 
formed. For example, due to the neutrinos in the final 
state, these decays are not ruled out by the limits on 
the branching ratios of lepton flavor- violating (LFV) de- 
cays t+ — > ^ ± 7T° at the 10 -8 — 10~ 7 level [35j since these 
searches require that the invariant mass of the £tt° system 
be consistent with m T . A loose requirement that the in- 
variant mass of the lepton-impostor pair be greater than 
1.5 GeV suppresses these branching ratios by a factor of 
about 4 x 10 -4 in the pseudoscalar case and 5 x 10~ 5 
in the scalar case, bringing them below the experimental 
limits on LFV decays. 

A difficulty in searching for these decays, especially for 
final states involving muons, is possible contamination 
from normal hadronic t+ — > p^v decays where p+ — > 
7r ± 7r° if the charged pion fakes a lepton. 



2. Associated Production 



> 

o 




3 4 5 6 

E sma\\ +E (*?,P,S) (GeV) 

FIG. 6. The cross sections for t + t~S (solid, red) with h T = 
1 and t + t~(tt°,P) (dashed, green) with = 1 and g T = 
1, with both t's decaying to electrons, and from t + t~ — ► 
e ± 7r T 7r c Vi/i' (dotted, blue) with n misidentified as e 2% of 



the time, shown as functions of E s 



+ E, 



The cut 



-^large ~\~ 

^(ir°,p,s) > -Ecm/2 = 5.29 GeV has been applied. 
Ismail is the smaller value of the charged lepton (or pion) 



energies and -Eiarge the larger. E, 



is the impostor or 



7T° energy. The r + r background has a cutoff at _Bsma.ii + 
E (7r°,P,s) = Ecm/2 = 5.29 GeV while the r+r" (tt , P, S) 
signals extend beyond that value. 

ton. To reduce such backgrounds, one can require that 

Ecm 



Plargc + E^a.p.s > 
-E-small + p7r°,P,S > 



2 ' 
Ecm 



(48) 
(49) 



One could also search for r pair production in asso- 
ciation with an impostor in e+e~ collisions. At y/s = 
10.58 GeV the cross sections for this are 

<r (e+e" t+t-tt ) ~ 9.8 {g*f pb, (45) 
a (e+e" -> t+t^P) ~ 9.8 g 2 T pb, (46) 
cr (e+e" t+t'S) ~ 110 h\ pb. (47) 

While these cross sections are large, it is difficult to sepa- 
rate the impostors or hardcore pion from the neutral pi- 
ons produced in hadronic tau decays. As far as we know, 
such searches have not been performed. The kinematics 
of these final states can offer some advantages in search- 
ing for impostors or a hardcore pion when compared with 
searches in r decays. 

A large background to t+t~ (tt , P, S) with both t's 
decaying to leptons comes from t+t~ production where 
one of the r's decays leptonically and the other to n^n® 
and a neutrino with the 7r+ being misidentified as a lep- 



where Pi ar ge (Psmaii) is the larger (smaller) of the charged 
lepton or pion energies. This cut ensures that the lep- 
tons (or charged pion faking a lepton) and the impos- 
tor (or hardcore pion) cannot have all originated from 
r decays following t+t~ production. To illustrate this, 
we generate t+t~ (tt , P, S) events with CalcHEP at 
\fs = 10.58 GeV with each r decaying to electrons, 
t+ — > e+Pi/, to minimize the pion misidentification rate. 
In order to demonstrate the effect of background from 
e+e~ — > t+t~ events, we also generate r+r~ events us- 
ing PYTHIA 6.4 [36| with one r decaying to an elec- 
tron, t — > e vi>, and the other to 7T+7T through a p, 
t —> p ± u — > ^+77°;/ and assign an electron misidentifi- 
cation rate of 2% due to charged pions. We require that 
-Eiarge + Pp.s > Pcm/2 and plot the resulting cross sec- 
tions as functions of P sma ii + Pp,s in Fig. [51 About 5% 
of t+t~P events satisfy both Eqs. ([48]) and {49]) as well 
as around 1% of t+t~ S events while none of the events 
from the r pair production background satisfy both of 
these cuts. After these cuts, the effective cross sections 
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for impostor or hardcore pion production become 



-T-7T 



P -> 

_ 5-> 



-4I/P) I 



4i/S)| 



15 

15 g 2 T fb, 
35 /i? fb. 



fb, 



(50) 
(51) 
(52) 



Using the best-fit couplings from the ir° transition form 
factor fits gives effective cross sections after these cuts 
in the range 0.5 — 25 fb, which would of course trans- 
late to O(10 3 ) events in the full BABAR dataset. Probing 
cross sections at this level will decisively test whether 
new (pseudo)scalars coupled to r are responsible for the 
observed rise in the ir° form factor. 

There is also a potential background from e + e~ — > 
77r + 7T~7r° with the photon traveling along the beam di- 
rection, going undetected. However, this can be elimi- 
nated by requiring missing transverse momentum. This 
also reduces many of the light quark-induced back- 
grounds generally. 

Charm pair production backgrounds, e.g., D + D~ pro- 
duction, are reduced by the cuts in Eqs. (|4"51) and (|4T)|) . 
just as the t + t~ background is. A potentially irreducible 
background comes from ir°D+D~ and ^D^D produc- 
tion with the D mesons decaying to leptons and neu- 
trinos, neutral pions, or to misidentified charged pions. 
We estimate the production cross section for these final 
states to be of the order 10 pb. The branching ratios 
to these final states are quite suppressed, especially to 
e^, which is another motivation for choosing t decays 
to electrons in this search. There is further suppression 
from the misidentification rate as well from the cuts in 
Eqs. (|4"5|) and (|4*g|) . causing these backgrounds to con- 
tribute at a level well below the relevant values of the 
signal cross section in Eqs. (|50)) - (|52|) . 



3. Effect on Hadronic r Decays 

The coupling of an impostor to t could lead to an 
apparent increase in the rate for — > ifiix^v, an ex- 
tremely important channel in estimating the hadronic 
spectral function used to calculate (g — 2) M . This could 
happen if a r emits an impostor in the three-body de- 
cay, r ± — > (P, S')7r ± i/. It could also occur if a r pair is 
produced in association with an impostor and one of the 
r's decays to ir^v. Taking both effects into account, we 
estimate a shift in this rate of 



ST (i 



r(r 



-> 7T°7r ± I/) 

r(r± ->• Tr ?!-^) 



0.26% g 2 



3.2% h 2 T . 



(53) 



(54) 



The shift in the hardcore pion case is more difficult to 
estimate, due to the interference with the standard decay 
amplitude. 

It is intriguing that this shift is not far from the level 
necessary to explain the longstanding discrepancy be- 
tween the hadronic spectral functions extracted from r 



decay and e + e~ data, see, e.g., |32j. This discrepancy 
is troubling since the hadronic spectral function pro- 
vides the largest theoretical uncertainty in the estimate 
of (g — 2) M . We leave a more complete analysis of this 
effect and the potential for it to alleviate the tension be- 
tween the different determinations of {g — 2) M for future 
work, but mention that such a solution appears possible. 

In addition, an increase in the di-pion branching of the 
t could change the measurement of the V us entry in the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix using the 
ratio of strange to non-strange r decays. Currently, the 
value of V us extracted in this way is 3.3cr to 3.6a smaller 
than the value predicted by the unitarity of the CKM 
matrix (341 . This sign of the change to V us looks 

correct in the models presented but more detailed study 
is required for quantitative conclusions. 



VI. CONCLUSIONS 

The BABAR measurement of the pion transition form 
factor did not agree with any of the QCD-based predic- 
tions that existed prior to the experimental publication. 
Experts' opinions are split whether the observed behav- 
ior can be accommodated within QCD even a posteri- 
ori, by, e.g., drastically modifying the pion distribution 
amplitude. Alternatively, quark models could provide a 
reasonable fit to data, but would require that the quark 
mass scale is on the order of 130 MeV, which is well above 
the current quark mass values. We have shown that this 
quark mass scale may be emergent as a nonperturbative 
correction due to the vacuum condensate-induced mix- 
ing between the massless quark and pion propagators, 
m° ff ~ — {qq)/(16F%). However, in order to explain the 
experimental rise of Q 2 J~(Q 2 ), one has to take the chiral 
quark model well outside its usual range of validity. Both 
explanations of the BABAR asymptotics of the pion form 
factor-QCD with a flattened distribution amplitude or 
a massive quark model-imply that the ir° field is more 
elementary than originally thought. 

The main focus of this paper is the new physics ex- 
planation of the form factor behavior, which is perhaps 
the least likely option among different logical possibili- 
ties outlined in the introduction. Nevertheless, we were 
able to construct some models of new physics that can be 
consistent with the observed form factor behavior. These 
models operate with new scalar or pseudoscalar particles 
coupled to t leptons, affecting T^q at the one-loop level. 
The choice of coupling to the r sector is motivated by 
necessity: e.g., similar models with light particles cou- 
pled to other leptons or to charm or bottom quarks are 
too constrained to play any role in the pion form factor. 
New particles can either imitate tt°, hence the name im- 
postors, or lead to the enhanced coupling of ir° to r's, 
thus supplying pions with a "hard core." Fitting exper- 
imental data, we were able to show that the range of 
couplings for such new states is consistent with exist- 
ing constraints, while the mass of the impostors must 



13 



be within a few MeV of m n . This might look like fine- 
tuning, except that similar accidents happen with SM 
particles (such as, e.g., the 5% coincidence of r lepton 
and D meson masses). The closeness of mi mpost or to the 
pion mass would naturally explain why there is no asso- 
ciated enhancement of the form factor for r\ and 77'. 

A relatively large coupling of new states to r leptons 
means that they can be searched for in a sizable r dataset. 
We have shown that both r decays, and especially asso- 
ciated production e + e~ — > t + t~(P, 5, it ) can be used 
to search for new "pion-like" physics. Choosing the op- 
timal couplings from the fits to the pion form factor, we 
find that B-factory samples should have on the order 
of O(10 3 ) events in the associated production channels 
with distinct kinematics that will facilitate their search. 
Such a large number of expected events gives hope that 
future experimental analyses can be conclusive: they ei- 
ther discover such new states or close perhaps the only 
new physics explanation of the unusual behavior of the 
pion transition form factor. 

Lastly, we note that interactions like those that we con- 
sidered could provide an explanation of the tension be- 



tween the hadronic vacuum polarization extracted from 
t and e + e~ data, which have important implications for 
the muon g — 2. Furthermore, there could be implications 
for the extraction of V us from r decays. Because of the 
importance of these subjects and detailed experimental 
data in the di-pion channels, further studies of hadronic 
r decays within these models might be warranted. 
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